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Summary  of  Research  Effort 


The  research  efforts  supported  by  this  grant  have  two  separate  but  closely 
related  thrusts:  (i)  Control  and  (ii)  parameter  estimation  (inverse  problems) 
for  distributed  parameter  systems  including  delay  and  partial  differential 
equations.  A  diverse  number  of  topics  have  been  investigated  and  can  be 
summarized  as  follows. 

I.  Control  (the  main  emphasis  was  the  development,  analysis,  and  numerical 
testing  of  computational  techniques): 

(i)  Numerical  methods  for  computation  of  feedback  gains  in  high  dimen¬ 
sional  LQR  problems  (a  hybrid  method:  Chandrasekhar  start  up;  modified 
Kleinman-Newton-Smith  iteration  for  the  gains). 

(ii)  Approximation  techniques  for  finite  dimensional  compensators  in  in¬ 
finite  dimensional  feedback  control  problems. 

(iii)  Approximate  solutions  of  operator  Riccati  equations:  regularity,  con¬ 
vergence. 

(iv)  Feedback  control  techniques  for  infinite  dimensional  systems. 

(v)  Control  and  stabilization  of  visco-elastic  structures. 

(vi)  Approximation  in  delay  and  Volterra  type  equations. 
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II.  Parameter  estimation  and  inverse  problems  (with  emphasis  on  develop¬ 
ment  and  analysis  of  computational  techniques  for  specific  applications  in 
mechanics  and  biology): 

(i)  Theoretical  frameworks  for  stability  and  convergence  of  approximation 
schemes:  variational  and  semigroup  formulations. 

(ii)  Regularization  techniques:  Compactness,  discretization,  Tikhonov 
techniques. 

(iii)  Augmented  Lagrangian  methods. 

(iv)  Estimation  of  damping  in  flexible  structures:  composite  and  vis¬ 
coelastic  materials. 

(v)  Estimation  of  growth,  predation,  dispersal,  size-structure  in  popula¬ 
tion  models. 

(vi)  Estimation  of  nonlinearites  in  distributed  systems. 
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On  a  Variational  Approach  to  Some  Parameter  Estimation  Problems 


We  consider  examples  (1-D  seismic,  large  flexible  structures,  bioturbation, 
nonlinear  population  dispersal)  in  which  a  variational  setting  can  provide  a 
convenient  framework  for  convergence  and  stability  arguments  in  parameter 
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A  Comparison  of  Stability  and  Convergence  Properties  of  Techniques  for 
.  Inverse  Problems  11 

Abstract 

We  consider  a  series  of  numerical  examples  and  compare  several  algorithms  for  estimation  of 
coefficients  in  differential  equation  models.  Unconstrained,  constrained  and  Tikhonov  regulariza¬ 
tion  methods  are  tested  for  the  behavior  with  regard  to  both  convergence  (of  approximation  meth¬ 
ods  for  the  states  and  parameters)  and  stability  (continuity  of  Lhe  estimates  obtained  with  respect 
to  perturbations  in  the  data  or  observed  states). 
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ABSTRACT 

We  discuss  techniques  for  the  estimation  of  nonlinearities  and  state-dependent 
coefficients  in  parabolic  partial  differential  equations.  Applications  to 
density-dependent  population  dispersal  and  nonlinear  growth/predation  models  are 
presented.  Computational  results  using  parallel  and  vector  architectures  are  discussed. 
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ABSTRACT 

We  discuss  techniques  for  the  estimation  of  nonlinearities  and  state-dependent 
coefficients  in  parabolic  partial  differential  equations.  Applications  to 
density-dependent  population  dispersal  and  nonlinear  growth/predation  models  are 
presented.  Computational  results  using  parallel  and  vector  architectures  are  discussed. 
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ESTIMATION  Of  STIFFNESS  AND  oaoMNC  IN  CANT  ILK  VEKtl)  KULCM-BKKNOULLI  BEAMS  WITH  TIP  BOD  IKS 
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Af>st  r.ict .  Wv  develop  finite  dimensional  approximation  schemes  for  tiie  Identification  of  spatially 
varying  material  parameters,  t.e.  flexural  stiffness  and  viscous  damping  coefficient!  in  hybrid  models 
for  flexible  beams  with  tip  bodies.  Our  schemes  are  derived  via  an  application  of  sp 1 1 ne- based  C.alerkln 
techniques  to  the  const* r vat i ve  form  state  space  representation  for  the  coupled  system  of  ordinary  ami 
partial  differential  equations  and  boundary  conditions  which  describe  tire  dynamics  of  the  system.  A 
convergence  theory  is  briefly  outlined  and  a  discussion  of  our  findings  based  upon  extensive  numerical 
studies  carried  out  on  Ih>iIi  conventional  and  vector  processors  is  included. 


Keywords.  Parameter  estimation;  flexible  structures;  damping;  approximation  theory;  computational 
methods;  splines. 


Ue  report  here  on  a  part  of  our 
Investigations  on  methods  to  estimate  continuum 
material  parameters,  and  specifically  stiffness 
and  damping  parameters  in  structures  described  by 
an  Euler  -  Bernoulli  formulation  of  the 
dynamics.  Our  efforts  have  been  motivated  in 
large  part  by  research  on  largo  f  texlhle  space 
structures.  Among  such  space  structures  of 
interest  are  plate  and  beam* like  truss 
structures,  space  platforms  with  multiple 
appendages  including  solar  panels,  robot  arms  and 
masts  with  tip  bodies  (e.g.  ilk*  shuttle  orblter 
with  beam/tlp  mass  conf 1 gurat ions )  and  large 
antennas.  These  structures  are  highly  flexible, 
possess  relatively  low  damping  characteristics, 
and  in  general  won't  support  their  own  weight  in 
earth's  gravity.  Their  flexible  members  (truss 
beams,  panels)  are  usually  constructed  of 
composite  materials  (e.g.  graph! Le  epoxy)  while 
membrane* 1 ike  woven  mesh  surfaces  are  found  in 
the  antenna  reflectors.  Thus,  in  these 
investigations  one  must  deal  with  complex 
composite  structures  with  variable  (in  apace  and 
time)  geometry  and  material  prone rc lex  for  which 
continuum  models  with  variable  parameters  offer 
•om:  obvious  advantages. 

The  Identification  or  paromeier  estimation 
techniques  which  we  describe  here  are  important 
for  model  development  and  analysis,  periodic 
auterlal  parameter  updates  in  structures,  as  well 
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as  being  a  precursor  to  end  Integral  part  of  the 
development  of  sophls C tested  control  and 
stabl l i tat  ion  strategies  for  large  flexible 
structures.  Ue  focus  here  on  the  possibilities 
of  estimating  viscoelastic  damping  of  Voigt* 
Kelvin  type  (Clough  and  Penzlen,  1975;  Popov, 
1968)  -  i.e.  where  the  damping  aomenl  is 
proportional  to  the  strain  rate.  While  there  is 
growing  evidence  that  damping  In  composite 
material  structures  Is  significantly  more  complex 
then  that  arising  from  the  Volgt-Kel vin 
hypothesis,  there  are  materials  for  which  the 
Voigt-Kelvln  assumption  is  a  good  approximation. 
Therefore  such  models  provide  an  excellent  class 
on  which  to  start  in  development  and  testing  of 
methods. 

In  addition  to  our  investigations  on 
approximation  methods  for  estimation  problems,  we 
have  also  been  exploring  the  use  of  special 
archl tectures-array  processors  snd  vector 
machines  *  in  conjunction  with  Che  methods  and 
algorithms  we  are  developing.  A  brief  discussion 
of  cenaln  aspects  of  our  findings  In  this  area 
Is  Included  below. 

The  results  presented  here  are  based  on  a 
conservative  state  apace  form  of  the  dynamical 
equations  for  a  cent  1  levered  beam  with  tip  body 
attached  to  the  free  end.  That  ^.s,  we  consider 
Che  higher  order  analogue  of  the  classical 
formulation  in  which  a  second  order  (in  time  and 
space  derivatives)  hyperbolic  equation  Is 
rewritten  as  a  first  order  vector  system  of 
equations  by  choosing  the  natural  states  of 
strain  ux  and  velocity  u(.  Due  to, space 
limitations,  we  shall  attempt  to  keep  details, 
discussions,  and  notation  to  a  minimum.  For  more 
complete  and  careful  discussions  we  refer  readers 
to  (banks  and  Rosen,  1985  and  1986). 

As  we  have  Indicated  above,  the  structures 
m  of  Interest  are  In  general  quite  complex.  The 
simplest  problems  Involve  beam  and  plate-like 
•odels  or  models  for  membrane- 1  Ike  surfaces*  One 
popular  and  useful  class  of  dynamical  system 
models  Involves  a  cantilevered  beam  with  tip  body 
and  base  acce lerat ion.  Such  models  are  useful, 
for  example,  in  studying  shuttle  deployed  masts 
with  attached  payload  -  e.g.  antenna  -  as 
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Computational  Methods  for  the  Identification  of  Spatially 
Varying  Stiffness  and  Damping  in  Beams 


ABSTRACT 


A  numerical  approximation  scheme  for  the  estimation  of  functional  parameters  in 
Euler-Bemoulli  models  for  the  tansverse  vibration  of  flexible  beams  with  tip  bodies  is  developed. 
The  method  permits  the  identification  of  spatially  varying  flexural  stiffness  and  Voigt-Kelvin 
viscoelastic  damping  coefficients  which  appear  in  the  hybrid  system  of  ordinary  and  partial 
differential  equations  and  boundary  conditions  describing  the  dynamics  of  such  structures.  An 
inverse  problem  is  formulated  as  a  least  squares  fit  to  data  subject  to  constraints  in  the  form  of  a 
vector  system  of  abstract  first  order  evolution  equations.  Spline-based  finite  element 
approximations  are  used  to  finite  dimensionalize  the  problem.  Theoretical  convergence  results  are 
given  and  numerical  studies  earned  out  on  both  conventional  (serial)  and  vector  computers  are 
discussed. 


THE  IDENTIFICATION  OF  A  DISTRIBUTED 
PARAMETER  MODEL  FOR  A  FLEXIBLE  STRUCTURE 


by 


H.T.  Banks,  S.S.  Gates,  I.G.  Rosen,  Y.  Wang 
August  1986  LCDS  #86-32 


ABSTRACT 


We  develop  a  computational  method  for  the  estimation  of  parame¬ 
ters  in  a  distributed  model  for  a  -flexible  structure.  The  structure 
we  consider  (part  of  the  " RPL  experiment")  consists  of  a  cantilevered 
beam  with  a  thruster  and  linear  accelerometer  at  the  free  end.  The 
thruster  is  fed  by  a  pressurized  hose  whose  horizontal  motion  effects 
the  transverse  vibration  of  the  beam.  We  use  the  Euler-Bernoulli 
theory  to  model  the  vibration  of  the  beam  and  treat  the  hose-thruster 
assembly  as  a  lumped  or  point  mass-dashpot-spring  system  at  the  tip. 
Using  measurements  of  linear  acceleration  at  the  tip,  we  estimate  the 
hose  parameters  (mass,  stiffness,  damping)  and  a  Voigt-Kelvin 
viscoelastic  structural  damping  parameter  for  the  beam  using  a  least 
squares  fit  to  the  data. 

We  consider  spline  based  approximations  to  the  hybrid  (coupled 
ordinary  and  partial  differential  equations)  system;  theoretical 
convergence  results  and  numerical  studies  with  both  simulation  and 
actual  experimental  data  obtained  from  the  structure  are  presented  and 
discussed. 
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A  Numerical  Algorithm  for  Optimal  Feedback  Gains 
in  High  Dimensional  LQR  Problems 


H.T.  Banks 
K.lto 


ABSTRACT 

We  propose  a  hybrid  method  for  computing  the  feedback  gains  in  linear 
quadratic  regulator  problems.  The  method,  which  combines  use  of  a  Chandrasekhar 
type  system  with  an  iteration  of  the  Newton-Kleinman  form  with  variable 
acceleration  parameter  Smith  schemes,  is  formulated  so  as  to  efficiently  compute 
directly  the  feedback  gains  rather  than  solutions  of  an  associated  Riccati  equation. 
The  hybrid  method  is  particularly  appropriate  when  used  with  large  dimensional 
systems  such  as  those  arising  in  approximating  infinite  dimensional  (distributed 
parameter)  control  systems  (c.g.,  those  governed  by  delay-differential  and  partial 
differential  equations).  Computational  advantages  of  our  proposed  algorithm  over  the 
standard  eigenvector  (Potter,  Laub-Schur)  based  techniques  are  discussed  and 
numerical  evidence  of  the  efficacy  of  our  ideas  presented. 


Kcv  Words:  LQR  problems,  feedback  gains,  distributed  parameter  systems, 
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proceeds  in  s  heuristic  and  iterative  fashion,  much  like  an  experimental 
inquiry  (see  Zia  1966  for  numerical  test  examples  and  details  on 
procedures  to  follow). 
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Analyzing  Field  Studies  of  Insect  Dispersal  Using  Two-Dimensional 
Transport  Equations 

ABSTRACT 

A  variety  of  diffusion  and  convection/diffusion  models  were  fit  to 
field  mark-recapture  data  (from  Havkes  1974)  for  female  cabbage  root 
flies  (Delia  brassicae  Bouche).  By  comparing  the  performance  of 
different  models  we  were  able  to  test  Havkes1  hypothesis  that  anemotaxis 
is  key  to  Delia*6  discovery  of  Brassies  crops.  Whereas  models  lacking  a 
convection  term  totally  failed,  models  with  a  convection  term  explained 
39-442  of  the  observed  variance  in  Del  is  recapture  density.  The 
direction  of  the  best-fit  convection  vector  va6  towards  Brassica  and  in 
almost  perfect  opposition  to  the  prevailing  winds.  This  suggests  that 
Dfel is  fly  upwind  in  the  presence  of  Brassica  odors.  The  application  of 
diffusion/ convection  models  to  insect  dispersal  is  discussed  in  general, 
with  special  emphasis  on  parameter  identification  methods.  These 
methods  (which  are  new  to  the  biological  literature)  allow  one  to  find 
realistic  dispersal  models  that  describe  field  mark-recapture  data. 
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Inverse  Problems  in  the  Modeling  of  Vibrations  of  Flexible  Beams 


ABSTRACT 

The  formulation  and  solution  of  inverse  problems  for  the  estimation  of  parameters  which 
describe  damping  and  other  dynamic  properties  in  distributed  models  for  the  vibration  of  flexible 
structures  is  considered.  Motivated  by  a  slewing  beam  experiment,  the  identification  of  a  nonlinear 
velocity  dependent  term  which  models  air  drag  damping  in  the  Euler-Bernoulli  equation  is 
investigated.  Galerkin  techniques  are  used  to  generate  finite  dimensional  approximations. 
Convergence  estimates  and  numerical  results  are  given.  The  modeling  of,  and  related  inverse 
problems  for  the  dynamics  of  a  high  pressure  hose  line  feeding  a  gas  thruster  actuator  at  the  tip  of  a 
cantilevered  beam  are  then  considered.  Approximation  and  convergence  are  discussed  and  numerical 
results  involving  experimental  data  are  presented. 
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Modeling  and  Estimation  in  Size  Structured  Population  Models 

ABSTRACT 

We  report  on  current  investigations  involving  both  deterministic  and  stochastic 
versions  of  size  structured  population  models  of  the  McKendrick  -  Von  Foerster 
type.  Simulation  studies  that  demonstrate  several  mechanisms  of  size  dispersion  are 
presented.  We  also  present  preliminary  computational  results  for  inverse  problems 
involving  the  estimation  of  survival  and  growth  parameters  for  larval  fish  from  size 
structured  field  data. 
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one  can  show  that  A(q)  generates  a  Cq  semigroup  T(t;q)  on  H.  If  b  >  0 
in  (B'),  this  semigroup  is  uniformly  exponentially  stable  and  if  in  (B')  we 
replace  U  by  I41  |  and  have  b  >  0,  the  semigroup  is  analytic.  For 
Euler— Bernoulli  beams,  the  general  case  handles  viscous  and  spatial  hysteresis 
damping  (with  uniform  exp.  stability  if  the  damping  coefficient  is  strictly 
positive)  while  Kelvlrr-Voigt  damping  is  included  in  the  analytic  semigroup 
case . 

Identification  problems  for  these  second  order  systems  may  be  formulated 
in  a  manner  analogous  to  the  first  order  case  outlined  above;  a  con¬ 
vergence/stability  theory  under  weak  compactness  assumptions  (typically  Q 
can  be  taken  as  a  subset  of  C(fi)  with  the  supremum  metric)  can  be  given 
using  the  resolvent  form  of  the  Trotter- Kato  theorem.  This  yields  results 
that  are  a  significant  improvement  over  those  currently  in  the  research 
literature.  Details  will  be  given  in  a  forthcoming  manuscript. 
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A  Theoretical  cramework  for  Convergence  and  Continuous  Dependence  of 
Estimates  in  Inverse  Problems  for  Distributed  Parameter  Systems 


Abstract 

In  this  note  we  announce  a  framework  in  which  one  can  treat  very 
general  classes  of  parameter  estimation  problems  for  distributed  systems.  Using 
this  approach  one  can  obtain  both  convergence  and  continuous  dependence 
(stability)  results  under  very  weak  regularity  and  compactness  assumptions  on 
the  set  of  admissible  parameters.  The  framework  includes  both  first  and  second 
order  abstract  systems  of  importance  in  applications. 
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Summary.  The  use  of  parameter  estimation  techniques  for 
partial  differential  equations  is  illustrated  using  a  predator- 
prey  model.  Whereas  ecologists  have  often  estimated  pa¬ 
rameters  in  models,  they  have  not  previously  been  able  to 
do  so  for  models  that  describe  interactions  in  heterogeneous 
environments.  The  techniques  we  describe  for  partial  differ¬ 
entia!  equations  will  be  generally  useful  for  models  of  inter¬ 
acting  species  in  spatially  complex  environments  and  for 
models  that  include  the  movement  of  organisms.  We  dem¬ 
onstrate  our  methods  using  field  data  from  a  ladybird  beetle 
( Coccinella  sepiempunciaia )  and  aphid  ( Urolcucon  nigroiu- 
herculatum)  interaction.  Our  parameter  estimation  algo¬ 
rithms  can  be  employed  to  identify  models  that  explain  bet¬ 
ter  than  80%  of  the  observed  variance  in  aphid  and  ladybird 
densities.  Such  parameter  estimation  techniques  can  bridge 
the  gap  between  detail-rich  experimental  studies  and  ab¬ 
stract  mathematical  models.  By  relating  the  particular  best- 
f>:  models  identified  from  our  experimental  data  to  Other 
information  on  Coccinella  behavior,  we  conclude  that  a 
term  describing  local  taxis  of  ladybirds  towards  prey  (aphids 
in  this  case)  is  needed  in  the  model. 

Key  words:  Parameter  estimation  -  Species  interaction  - 
redistribution  -  Heterogeneous  environments 


Although  mathematical  models  of  species  interactions  com¬ 
mand  much  attention  in  the  ecological  literature,  their  use¬ 
fulness  is  sometimes  questioned  (see.  for  example.  SimberiofT 
19S3).  One  problem  with  models  in  ecology  is  that  they 
are  difficult  to  explicitly  relate  to  field  data.  By  assuming 
that  interacting  populations  are  at  an  equilibrium.  Schoener 
(1974)  and  Belovskv  (1984)  have  been  able  to  use  numerical 
least-squares  methods  to  estimate  parameters  in  models  on 
the  basis  of  field  data.  These  efforts,  however,  depend  on 
explicit  analytic  representations  for  the  functions  (solutions) 
being  observed  or  fit  to  data.  Usually  such  explicit  solutions 
arc  possible  only  for  the  simplest  of  models  (Mueller  and 
Ayala  1981)  or  for  the  equilibrium  points  of  more  complex 
models.  Furthermore,  in  the  ecological  literature  parameter 
estimation  techniques  have  not.  to  our  knowledge,  been  ap¬ 
plied  to  transient  dynamic  models  of  spatially  structured 
interactions  (i.e.,  interactions  in  which  densities  vary  in 
space  and  are  altered  by  dispersal).  Since  there  is  a  growing 


interest  in  non-equilibrium  models  (Chesson  and  Case  1985) 
and  in  models  that  allow  organisms  to  move  about  in  a 
heterogeneous  environment  (Levin  1974,  Wiens  et  al.  1985). 
there  is  a  need  for  parameter  estimation  techniques  that 
do  not  rely  on  writing  down  a  particular  formula  which 
represents  the  solution  to  a  model. 

In  this  paper,  we  present  parameter  estimation  tech¬ 
niques  for  models  that  take  the  form  of  partial  differential 
equations.  Such  models  can  be  used  to  represent  species 
interactions  tn  heterogenous  environments  plus  a  wide  vari¬ 
ety  of  movement  behaviors.  The  mathematical  theory  be¬ 
hind  these  techniques  is  developed  elsewhere  (Banks  and 
Kareiva  1983;  Banks  et  al.  1985;  Banks  and  Murphy  1986; 
1987).  Our  goal  here  is  to  make  experimental  ecologists 
aware  that  if  they  collect  spatially  and  temporally  structured 
field  data  (i.e.,  densities  at  position  x,.  x2....,  xm  at  each 
of  the  times  r„),  there  is  the  possibility  of  using 

these  data  to  identify  parameters  in  models  of  interacting 
populations  with  dispersal.  We  will  illustrate  the  use  of  such 
parameter  estimation  techniques  by  applying  our  methods 
to  field  data  we  have  collected  on  a  predator-prey  interac¬ 
tion  between  ladybird  beetles  (Coccinella  sepiempunciaia) 
anc  aphids  (Urolcucon  nigrotuberculaium). 


The  parameter  estimation  problem 

for  general  interaction  and  redistribution  models 

Interaction  and  redistribution  of  species  can  be  generally 
modeled  by  equations  of  the  form 
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where  A’)  anc  A’;  are  the  densities  of  two  interacting  species 
that  vary  in  space  (i.e..  with  respect  to  x).  The  /,  and  /, 
represent  the  local  "h.nctics''  or  "interaction  terms'  ior 
each  species  -  these  could  be  competition  functions,  preda¬ 
tor-prey  terms,  or  host -disease  terms.  Tne  first  two  terms 
on  the  right  side  of  fla)  and  (lb)  describe  the  movement 
(or  "redistribution”)  of  A',  and  A';  with  respect  to  x.  The 
movement  terms  with  D,  in  the  parentheses  correspond  to 
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eters  for:  nonlinear  interaction  terms,  nonlinear  diffusion 
terms,  nonlinear  convection  terms,  and  taxis  towards  or 
away  from  other  organisms.  Thus  we  will  clearly  have  the 
ability  to  identify  biologically  realistic  models  of  species  in¬ 
teractions  in  spatially  heterogeneous  environments. 

Our  approach  using  parameter  estimation  is  not.  how¬ 
ever.  a  panacea  for  population  modelers.  Our  experience 
indicates  that  the  models  which  describe  field  data  well  usu¬ 
ally  include  numerous  parameters  and  are  difficult  to  ana¬ 
lyze  (Banks  and  Kareiva  1983:  Banks  ct  al.  I985|.  Some¬ 
times  we  are  left  w-ith  detailed  models  that  can  only  be 
explored  by  exhaustive  numerical  studies.  In  addition,  there 
is  the  danger  that  so-called  "best-fit''  models  are  biological 
nonsense.  The  techniques  we  have  used  are  so  powerful 
that  rarely  do  we  fail  to  find  some  model  which  generates 
numerical  solutions  closely  matched  to  the  data. 

Parameter  and  model  identification  approaches  are  best 
viewed  as  one  tool  in  the  arsenal  of  population  biologists. 
Once  population  data  have  led  to  a  suite  of  models  such 
as  examples  A-D  in  Table  1.  experiments  need  to  be  per¬ 
formed  to  test  the  models  or  components  of  the  models. 
For  instance,  using  the  equations  presented  in  Table  1  we 
could  make  predictions  about  ladybird  and  aphid  popula¬ 
tion  dynamics  for  any  set  of  initial  conditions,  and  then 
lest  those  predictions.  Alternatively,  theoreticians,  more  in¬ 
terested  in  general  features  of  models,  could  examine  the 
properties  of  models  that  are  simplified,  but  qualitatively 
similar  to  those  identified  in  Table  1.  , 

Ideally,  we  would  like  to  verify  models  by  direct  observa¬ 
tion  of  individual  demography  and  individual  movement 
behavior.  Unfortunately,  these  sorts  of  observations  are  not 
always  possible.  Indeed,  we  were  able  to  make  these  obser¬ 
vations  on  Coccinella  only  when  we  arranged  a  0.1  m  x  10  m 
line  of  goldenrod  as  our  experimental  arena.  As  the  vegeta¬ 
tion  becomes  more  two-dimensional  (even  a  one-meter  wide 
strip),  detailed  behavioral  observations  become  increasingly 
difficult.  In  a  wide  variety  of  settings  (e.g.  pest  management, 
studies  of  the  spread  of  invading  organisms,  etc.),  data  on 
population  densities  at  different  times  and  places  may  be 
the  only  information  that  is  practically  available.  In  these 
cases,  parameter  estimation  algorithms  can  help  generate 
reasonable  first-guesses  for  models  of  population  interac¬ 
tion,  which  can  then  focus  subsequent  empirical  efTorts. 
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ABSTRACT 

A  convergence  theory  is  given  for  approximation  techniques  to  treat  inverse  problems 
involving  systems  of  nonlinear  parabolic  partial  differential  equations.  These  techniques  can 
be  used  to  estimate  density-dependent  dispersal  coefficients  in  population  models,  as  well  as 
nonlinear  growth  and  predation  terms.  Numerical  experiences  with  the  resulting  algorithms 
on  both  conventional  (scalar)  and  vector  computers  are  reported  along  with  an  indication  of 
performance  of  the  methods  with  field  data  from  prey-predator  experiments. 
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We  survey  some  problems  and  related  recent  results  in  computational  methods  for  two 
classes  ot  problems  for  distributed  parameter  systems-  parameter  estimation  or  inverse 
problems  and  feedback  control. 


I.  INTRODUCTION 


In  this  lecture  wc  shall  give  a  survey  of  some  of  our 
recent  and  current  efforts  in  the  Center  for  Control 
Sciences  at  Brown  University.  Wc  shall  restrict  our 
attention  to  research  projects  involving  theoretically 
sound  computational  aspects  ot  parameter 
estimation  and  feedback  control  in  problems 
described  by  partial  differential  equations.  As  lime 
permits,  the  topics  listed  below  will  be  discussed. 

II.  THEORETICAL  FRAMEWORK  FOR 

parameter  estimation 

Wc  consider  first  order  systems  .dependent  on 
parameters  qt Q  as  given  by  an  abstract  equation 

u(t)=A(q)u(t)+F(t,q) 

u(0)=uo(q) 

in  a  Hilbert  space  H.  We  scek_  to  estimate  the 
parameters  q  using  observations  u-,  for  u(t;;q)  in  a 
least-squares  fit  criterion.  We  assume  that  Q  is  a 
metric  space  with  metric  d  and  for  qrQ,  A(q)  is  the 
infinitesimal  generator  of  a  Cg  semigroup  on  H.  We 
assume  that  A(q)  is  defined  through  a  parameter 
dependent  sesquilincar  form  o(q)(y)  associated  with 
the  weak  or  variational  form  of  (2.1).  Briefly,  let  V 
be  a  Hilbert  space  that  is  continuously  and  densely 
imbedded  in  H.  Denote  a  family  of  parameter 
dependent  sesquilincar  forms  by  o(q):VxV-C,qrQ. 
Wc  assume  that  o  possesses  the  following 
properties: 

(A)  Continuity:  For  q,  q  «Q,  wc  have  for  all  q,  4*  V 
|o(q)((t>,^)- o(q)(q,^)!  <  <J(q.q>I«S>lv  I’f’lv- 


(B)  Coercivitv:  There  exists  ct  >0  and  some  i  such 
that  for  q«Q,C(  V  we  have 


(C)  Boundedness:  There  exist  Cj  >0  such  that  for 
qrQ,  f\<i«  V  wc  have 

|o(q)(q.iS)|<  Cj|e|v|»!<|v. 

Under  these  assumptions,  a  defines  in  the  usual 
manner  operators  A(q)  such  that  o(q)(0,«J/)  = 

<-A(q)C,<!<>  H  for  4>  c  dom  (A(q)),  4>  <  V  with 
dom(A(a))  dense  in  V.  Furthermore,  A(q)  is  the 
generator  of  an  analytic  semigroup  T(t;q)  on  H 
(indeed,  A(q)  is  sectorial  with  f  xl-A(q))dom(A(q)) 

=  H).  Property  (B)  guarantees  that  the  resolvent 
operator  Rv(A(q))  s^IA^q))"1  exists  as  a  bounded 
operator  on  H. 

We  may  then  consider  Galerkin  type  approximations 
in  the  context  of  these  sesquilincar  forms.  Let  HN  be 
a  family  of  finite  dimensional  subspaces  of  H 
satisfying  PN’z-z.  for  zcH  where  PN  is  the  orthogonal 
projection  of  H  onto  HK.  We  further  assume  that 
HkcV  and  possesses  certain  V-approximaiion 
properties.  It  wc  now  consider  the  restriction  of 
o(q)(y)  to  HN’,  we  obtain  operators  AK(q)-HN-HN 
which,  because  of  (B),  satisfy  a  uniform  dissipative 
inequality  and  can  be  shown  to  generate  semigroups 
T^itq)  in  HN'. 

Using  a  resolvent  convergence  form  of  the  Trotter- 
Kato°thcorcm  to  establish  convergence  properties  of 
TN(t;qN)  to  T(t;q)  whenever  {qR}  is  any  arbitrary 
sequence  converging  to  q  in  d-metric,  one  can 
establish  a  convergence  and  stability  theory'  for  the 
resulting  approximating  least-squares  estimation 
problems.  ' 

Among  the  examples  that  can  be  treated  immediately 
with  the  above  theory  are  the  usual  coefficient  and 
boundary  parameter  identification  problems  for 
parabolic  systems  (in  this  case  H  =  1^(0)  and 
V  =  Ho(n)  or  some  other  appropriately  chosen 
subspace  of  HJ(n)  with  modified  boundary 
conditions).  This  theoretical  framework  can  also 
readily  be  used  for  domain  identification  problems 
arising  in  thermal  tomography. 


The  first  step  of  our  hybrid  method  requires  the 
solution  of  n(m+p)  simultaneous  equations,  while 
each  iteration  of  the  usual  Newton-Kieinman  step 
requires  the  solution  of  a  Lyapunov  equation  for  the 
nxn  symmetric  estimates  of  P. 

However,  one  can  use  factorization  ideas  and  the 
Smith  method  for  Lyapunov  equations  to  reformulate 
the  Newton-KJeinman  method  as  a  direct  iterative 
method  for  the  mxn  gain  K,  thereby  providing 
additional  computational  advantages.  To  speed  up 
our  calculations  and  improve  convergence  in  the 
Smith  algorithm,  we  have  used  a  variable  stepsize 
Smith  method  to  solve  the  Lyapunov  equations. 

Our  initial  numerical  experiments  with  this  hybrid 
algorithm  have  provided  quite  encouraging  results. 
For  example,  if  n  is  the  dimension  of  our 
approximating  systems,  we  have  observed  in  our 
examples  with  1-dimensional  parabolic  systems  that 
the  traditional  eigenvector  based  methods  (Potter, 
Laub-Schur)  exhibit  computational  effort  (and  time) 
that  grows  like  n3  whereas  the  hybrid  method  time 
grows  like  n. 


IV.  CONCLUDING  REMARKS 

Much  of  the  research  discussed  in  this  survey  lecture 
involves  very  recent  and  current  results.  Some  of  it 
involves  joint  efforts  with  other  investigators 
(including  K.  I  to,  Y.  Wang,  F.  Kojima,  DJ.  Inman,  H. 
Cudney,  R.  Fabiano,  G-  Rosen,  S.  Reich  and  C 
Wang)  and  in  several  cases  manuscripts  with  the 
detailed  results  are  currently  being  written.  At  the 
time  of  the  lecture  we  shall  provide  a  bibliography 
available  on  these  and  related  results. 
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Abstract 

«C  report  on  our  use  of  spline  based  inverse 
procedures  to  estimate  damping  coefficients  in  distributed 
parameter  systems  for  fleaible  structures.  Damping  models 
involving  viscous  (air)  damping  and  Kelvin-Voigt  damping 
in  an  Eulct-Bcrnoulli  framework  are  used  to  analyse  data 
from  vibration  experiments  with  composite  material  beams. 

1.  Intrnduct ion 

In  this  presentation  we  report  on  continuing  efforts 
that  we  love  pursued  for  the  past  year.  These  efforts 
involve  a  combination  of  experimental  investigations  with 
the  simultaneous  development  of  the  mathematical  and 
computational  aspects  of  a  theoretical  methodology  to 
support  the  experiments. 

Our  long  term  goals  include  the  understanding  of 

damping  mechanisms  in  complex  distributed  structures 
constructed  from  composite  materials.  Our  ouest  has 

obvious  motivation  from  and  potential  consequences  for  the 
design  of  control  (active  and  passive)  systems  for  large 

flexible  space  structures. 

The  initial  efforts  on  which  we  report  here  email  the 
testing  ard  development  of  spline  based  approximation 
technioues  for  inverse  problems  arising  in  attempts  to 
model  and  quantify  damping  in  composite  material  beams. 
Substantial  previous  mathematical  efforts  along  with 
numerical  tests  on  simulated  data  can  be  found  in  the 

literature  (see  JBCRJ.  fBR)  for  a  sample  of  some  of  these 
results  and  references  to  other  work). 

We  briefly  outline  the  underlying  mathematical  ideas, 
describe  the  experimental  equipment  employed,  and  then 
report  on  two  representative  experiments  and  our  findings 
using  our  techniques  with  the  data  from  these  experiments. 

IL  Mathematical  Foundations 


The  basic  mathematical  model  that  we  have 
considered  in  connection  with  the  efforts  on  which  we 
report  in  this  note  consists  of  the  equations  for  a 
cantilevered  Euicr-BernouUi  equation  with  tip  mass, 
damping,  and  an  apDlied  transverse  force.  Specifically  we 
take  as  our  model  the  system  for  transverse  vibrations  of  a 
long  slender  (length  I)  cantilevered  beam: 


p  ?-¥-  ♦  El 
ex' 


2lu.  .  c  i 

ox*  c  cx'di 


*  7’ 


.011 


(2.1) 


■  f(t,xp.T),  0  <  x  <  f.  t  >  0, 


u(u0)  -  -“MU))  -  0,  t  >  0 
ox 


(^2) 


El  f“-  (L •)  ♦  CpI  "“(t.O  -  0.  I  >  0. 

nr *  v  nr*rT 


(2J) 


El  ((.l)  *  c— I  (t.l)  -  m  ((./),  t  >  0,  (Z4) 

ex'*  v  cx*e i  c:‘ 


u(0.x)  -  d(x).  (0,x)  -  *<x),  0  <  x  <  1.  CL5) 

Cl 


Here  p  is  the  linear  mass  density,  m  is  the  tip  mass,  cDl 
and  7  arc  damping  coefficients  for  Kelvin-Voigt  and 
viscous  (air)  damping,  respectively,  and  f  is  a  parameter 
dependent  (xp  (or  the  point  of  impact,  T  for  the  duration) 
term  representing  a  transverse  force  such  as  that  obtained 
with  a  force  hammer.  The  functions  t>  and  4>  represent 
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given  initial  displacement  and  velocity  of  the  beam. 

We  assume  that  we  arc  given  observations  -  in  this 
case,  of  lip  acceleration  -  from  which  wc  wish  to  estimate 
unknown  parameters.  In  general  these  include  the  stiffness 
coefficient  El  as  well  as  the  damping  parameters  c^l.  7 
and  the  input  parameter  T.  Denoting  this  collection  of 
unknown  parameters  by  Q  ■  (E1.cd1,7.t),  wc  may  formulate 
this  inverse  problem  as  the  leasi-SQuarcs  estimation  problem 
of  minimizing  the  functional 

J(q)  -  I  (i,.»:q)  -  »(«,.) I ’  (2.6) 

I  1  Ol  1 

over  an  admissible  class  Q  of  parameter  values.  In  this 
fii-to-data  criterion  (a(t,))  are  the  observations  (data)  for 
the  tip  acceleration  (at  a  »  J)  of  the  beam. 

This  estimation  problem  cannot,  of  course,  be  solved 
analytically:  however,  one  can  employ  iterative  optimization 
schemes  coupled  with  an  approximation  method  for  the 
infinite  dimensional  system  (2.1)  -  (1.5).  For  the  efforts 
discussed  here,  we  approximated  the  system  via  Galerktn 
procedures  using  cubic  B~sp)ine  elements  as  explained  in 
(BGRW).  This  results  in  an  approximate  solution  uK  which 
satisfies  a  readily-solved  finite  dimensional  system 
approximating  (2-1)  -  (L_5).  We  use  this  in  the  criterion 
(2.6)  in  place  of  u,  thereby  obtaining  a  family  of 
approximating  estimation  problems  involving  minimization 
over  Q  of 

jn(q)  -  r  a?) 

1  *  ol  1 

Solving  these,  w-  obtain  a  scoucncc  of  estimates  fq*)  f or 
best-fit  parameters  in  the  original  estimation  problem.  A 
theory  for  convergence  of  these  estimates  can  be  obtained 
(see  c-&~  IBCRyBR)JBGRW]  for  details),  but  we  shall  not 
discuss  this  aspect  of  the  computational  techniques  here 
since  the  focus  of  this  report  is  the  jug  of  these 
techniques  in  studying  damping  mechanisms  in  composite 
material  structures. 

Wc  have  assumed  here  that  the  observations  arc  taken 
with  tip  accelerometers  since  that  is  precisely  the  case  for 
the  experiments  discussed  below.  However,  wc  could 
equally  well  assume  velocity  observations  are  taken;  indeed, 
wc  have  also  used  these  methods  with  experiments  in 
which  a  laser  vibrometcr  was  used  to  obtain  velocity 
measurements.  These  and  related  findings  will  be  detailed 
in  a  separate  manuscript  currently  in  preparation. 

Furthermore,  in  our  findings  discussed  here,  we  have 
concentrated  on  attempts  to  characterize  the  damping  in 
terms  of  Kelvin-Voigt  and/or  viscous  (air)  damping.  There 
is  considerable  evidence  to  suggest  that  damping 
mechanisms  in  composite  materials  arc  significantly  more 
complex  than  the  ones  described  by  a  combination  of 
viscous  and  Kelvin-Voigt  damping.  Indeed,  current 
conjectures  include  models  that  involve  hysteretic  or 
hereditary  (time  and/or  space)  effects.  Recent  theoretical 
results  fBI]  provide  a  convenient  mathematical  framework 
for  treating  inverse  problems  for  such  models;  we  arc 
currently  using  the  associated  computational  techniques  to 
investigate  the  possibility  of  describing  experimental  data 
in  terms  of  models  with  hystcrctic  damping. 

lit.  Experimental  Procedures 

To  test  the  above  described  estimation  procedures  and 
initiate  our  damping  mechanism  studies  a  scries  of 
experiments  were  carried  out  at  the  Mechanical  Systems 
Laboratory  located  at  SUNYaE.  Cantilevered  composite 
beams  with  a  removable  tip  mass  were  used  as  the  test 
structures.  Each  beam  consisted  of  a  one  meter  long  piece 
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Figure  6 


If  we  try  to  estimate  ^  and  cDI  using  modal 
techniques  (via  %  critical  damping  and  principles  relating  ( 
to  >  and  cDl  similar  to  that  noted  in  the  discussion  of  the 
previous  experiment),  we  obtain  values  7  .  -  .2982  slug/(fi 
sec)  and  c0Ics)  -  .4511  slug  fts/sec.  If  we  use  these  in 
the  model  equations  (2.1)  -  (2.5).  we  find  that  the  response 
is  damped  out  much  too  quickly  and  does  not  resemble 
the  experimental  data. 

A  final  comment  on  our  findings  in  this  example  has 
important  implications  for  use  of  modal  methods  in 
attempts  to  understand  damping  mechanisms  in  distributed 
parameter  structures.  While  we  found  that  the  viscous 
damping  (7)  was  extremely  important  in  the  overall 
modeling  of  the  low  frequency  mode  and  the  Kelvin-Voigt 
damping  (cpI)  was  most  important  in  modeling  correctly 
the  two  high  frequency  modes,  the  damping  does  not 
decouple  on  the  modes.  We  bclivc  this  poses  serious 
difficulties  for  any  attempts  to  understand  damping 
mechanisms  using  only  modal  techniques. 
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Grant  AFOSR-85-01 19  (DJ.I.). 
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An  Approximation  Theory  for  the  Identification  of  Nonlinear 
Distributed  Parameter  Systems 

ABSTRACT 

An  abstract  approximation  framework  for  the  identification  of 
nonlinear  distributed  parameter  systems  is  developed.  Inverse  problems  for 
nonlinear  systems  governed  by  strongly  maximal  monotone  operators 
(satisfying  a  mild  continuous  dependence  condition  with  respect  to  the 
unknown  parameters  to  be  identified)  are  treated.  Convergence  of 
Galcrkin  approximations  and  the  corresponding  solutions  of  finite 
dimensional  approximating  identification  problems  to  a  solution  of  the 
original  infinite  dimensional  identification  problem  is  demonstrated  using 
the  theory  of  nonlinear  evolution  systems  and  a  nonlinear  analog  of  the 
Trottcr-Kato  approximation  result  for  semigroups  of  bounded  linear 
operators.  The  nonlinear  theory  developed  here  is  shown  to  subsume  an 
existing  linear  theory  as  a  special  case.  It  is  also  shown  to  be  applicable 
to  a  broad  class  of  nonlinear  elliptic  operators  and  the  corresponding 
nonlinear  parabolic  partial  differential  equations  to  which  they  lead.  An 
application  of  the  theory  to  a  quasilinear  model  for  heat  conduction  or 


mass  transfer  is  discussed. 
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Semigroup  Theory  and  Numerical  Approximation  for  Equations  in  Linear 
Viscoelasticity 


Abstract 


We  consider  the  following  abstract  integro-differential  equation 


u(0~-4 


Eu(t)  —  J  g(s)u{t  +  s)  ds 


=  /(  0 


on  a  Hilbert  space.  Such  equations  arise  in  the  modeling  of  linear  viscoelastic  beams. 
The  equation  is  reformulated  as  an  abstract  Cauchy  problem,  and  several  approximation 
schemes  are  discussed.  Well-posedness  and  convergence  results  are  given  in  the  context  of 
linear  semigroup  theory.  Results  of  numerical  eigenvalue  calculations  for  various  approxi¬ 
mation  schemes  are  discussed. 
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,  Augmented  Lagrangian  Method  for  Equality  and  Inequality  Constraints  in 
Hilbert  Spaces 

Abstract 

In  this  paper,  we  consider  a  class  of  Lagrange  multiplier  methods,  called  the 
augmented  Lagrangian  method  for  the  minimizations  with  equality  and  inequality 
constraints  in  Hilbert  spaces.  We  obtain  a  local  square-root  convergence  of  this 
method  without  strict  complementary  conditions  when  finite  many  inequality 
constraints  are  augmented.  This  result  seems  to  even  generalize  the  known  results 
in  the  finite  dimensional  case. 


On  the  Regularity  of  Solutions  of  an  Operator  Riccati 
Equation  Arising  in  Linear  Quadratic  Optimal  Control 
Problems  for  Hcrediatry  Differential  Systems 

by 

Kaaufumi  lio 

May  1 9S6  LCDS 


OK  THE  REGULARITY  OF  SOLUTIONS  OF 


AN  OPERATOR  RICCATI  EQUATION  ARISING 
IN  LINEAR  QUADRATIC  OPTIMAL  CONTROL 
PROBLEMS  FOR  HEREDITARY  DIFFERENTIAL  SYSTEMS* 


by 


Karufumi  I to 


Lefschet:  Center  for  Dynamical  Systems 
Division  of  Applied  Mathematics 
Brown  University 
Providence,  RI  02912 


May  1986 


This  research  was  supported  by  the  Air  Force  Office  of  Scientif 
Research  under  grants  AFOSR -84-0398  and  AFOSR-E5-0305,  and  the 
National  Aeronautics  and  Space  Administration  under  grant  NAG-1 


Orr  the  Regularity  of  Solutions  of  an  Operator  Riccati  Equation  Arising  in 
Linear  Quadratic  Optimal  Control  Problems  for  Hereditary  Differential  Systems 


Abstract 

In  this  paper  operator  Riccati  equations  of  evolution  and 

algebriac  types  arising  in  the  linear  quadratic  optimal  control 

problem  for  hereditary  differential  systems  are  considered.  Regularity 

results  of  their  solutions  are  obtained.  Spectrally,  for  the  system 

with  single  point  delay  we  show  that  the  optimal  feedback  kernel 
10  ‘  00 

function  T.  is  C  .  This  study  is  motivated  by  our  numerical 
experiments  IS]  in  computing  the  optimal  feedback  gain  via  an 
algebraic  Riccati  equation,  which  indicate  that  such  a  function  is 
smooth.  The  regularity  result  is  essential  in  obtaining  a  rate  of 
convergence  of  numerical  approximations  of  the  optimal  feedback  gain 
operator. 
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"Approximation  of  infinite  delay  and  Volterra  type  equations" 
by  K.  Ito  and  F.  Kappel 


Summary 


Linear  autonomous  functional  differential  equations  of  neutral  type 
include  Volterra  integral  and  integrocif f erential  equations  as  special 
cases.  The  paper  considers  numerical  approximation  of  solutions  to 
these  equations  by  first  converting  the  initial  value  problem  to  an 
abstract  Cauchy  problem  in  a  product  space  0Rn  *  weighted  L^-space)  and 
then  using  abstract  approximation  results  for  C  -semigroups  combined  with 
Galerkin  type  ideas.  In  order  to  obtain  concrete  schemes  subspaces  of 
Legendre  and  Laguerre  polynomials  are  used.  The  convergence  properties 
of  the  algorithms  are  demonstrated  by  several  examples. 
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Finite  Dimensional  Compensators  for  Infinite  Dimensional  Systems  via 
Galerkin-Type  Approximation 


Abstract 

In  this  paper  we  discuss  existence  and  construction  of  stabilizing 
compensators  for  linear  time-invariant  system  defined  on  Hilbert  spaces.  We 
establish  an  existence  result  using  Galcrkin-type  approximation  in  which 
independent  basis  elements  are  used  instead  of  the  complete  set  of  eigenvectors. 
A  design  procedure  based  on  approximate  solutions  of  optimal  regulator  and 
optimal  observer  via  Galerkin-type  approximation  is  given  and  the  Schumacher 
approach  is  used  to  reduce  a  dimension  of  compensators.  A  detail  discussion 
for  parabolic  and  hereditary  differential  systems  is  included. 
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Abstract 

In  this  paper  a  new  technique  for  the  estimation  of  parameters  in  elliptic  partial 
differential  equations  is  developed.  It  is  a  hybrid  method  combining  the  output  *  least  - 
square  and  the  equation  error  method.  The  new  method  is  realized  by  an  augmented 
Lagrangian  formulation  and  convergence  as  well  as  rate  of  convergence  proofs  are  provided. 
Technically  the  critical  step  is  the  verification  of  a  coercivity  estimate  of  an  appropriately 
defined  Lagrangian  functional.  To  obtain  this  coercivity'  estimate  a  seminorm  regularization 
technique  is  used. 
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•  Chandrasekar  Equations  for  Infinite  Dimensional  Systems:  Part  II.  Unbounded 
Input  and  Output  Case 

Abstract 

A  set  of  equations  known  as  Chandrasekhar  equations  arising  in  the  linear 
quadratic  optimal  control  problem  is  considered.  In  this  paper,  we  consider  the 
linear  time-invariant  systems  defined  in  Hilbert  spaces  involving  unbounded  input  and 
output  operators.  For  a  general  class  of  such  systems,  we  derive  the  Chandrasekhar 
equations  and  establish  the  existence,  uniqueness,  and  regularity  results  of  their 


solutions. 
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Strong  Convergence  and  Convergence  Rates  of  Approximating  Solutions 
for  Algebraic  Riccati  Equations  in  Hilbert  Spaces 


6  bairasi 

In  this  paper,  we  consider  the  linear  quadratic  optimal  control 
problem  on  infinite  time  interval  for  linear  time-invariant  systems  defined  on 
Hilbert  spaces.  The  optimal  control  is  given  by  a  feedback  form  in  terms  of 
solution  n  to  the  associated  algebraic  Riccati  equation  (ARE).  A  Ritz  type 
approximation  is  used  to  obtain  a  sequence  nN  of  finite  dimensional 
approximations  of  the  solution  to  ARE.  A  sufficient  condition  that  shows  nN 
converges  strongly  to  n  is  obtained.  Under  this  condition,  we  derive  a 
formula  which  can  be  used  to  obtain  a  rate  of  convergence  of  nN'  to  a 
We  demonstrate  and  apply  the  results  for  the  Galcrkin  approximation  for 
parabolic  systems  and  the  averaging  approximation  for  hereditary  differential 


systems. 
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CHANDRASEKHAR  EQUATIONS  FOR  INFINITE  DIMENSIONAL 

SYSTEMS* 

KAZUFUM1  ITO*  and  ROBERT  K.  POWERS: 


Abstract.  In  this  paper  we  derive  the  Chandrasekhar  equations  for  linear  time  invariant  systems  defined 
on  Hilbert  spaces  using  a  functional  analytic  technique.  An  important  consequence  or  this  is  that  the  solution 
to  the  evolutional  Rjccati  equation  is  strongly  differentiable  in  time  and  one  can  define  a  "strong"  solution 
of  the  Rjccati  differential  equation.  A  detailed  discussion  on  the  linear  quadratic  optimal  problem  for 
hereditary  differential  systems  is  also  included. 


Key  words.  Chandrasekhar  equations,  Rjccati  operator,  regularity  results,  infinite  dimensional  systems 
AMSfMOS)  subject  classification.  49 

1.  Introduction.  The  Chandrasekhar  equations  [14]  are  an  alternative  form  to  the 
Rjccati  equations  from  which  the  optimal  feedback  gain  operator  may  be  calculated 
directly.  If  the  system  has  a  small  number  of  inputs  and  outputs,  the  Chandrasekhar 
algorithm  offers  significant  reduction  in  the  computational  complexity  for  determining 
the  optimal  feedback  gain.  As  observed  in  [20],  this  is  much  more  evident  in  the  infinite 
dimensional  case  if  the  optimal  feedback  gain  operator  is  calculated  numerically  using 
some  approximation  method.  In  this  case,  the  number  of  states  grows  linearly  to  the 
order  of  approximation. 

The  purpose  of  this  paper  is  to  derive  Chandrasekhar  equations  for  systems  defined 
by  evolution  equations  on  Hilbert  spaces  in  which  the  input  and  output  operators  are 
assumed  to  be  bounded.  The  form  of  the  Chandrasekhar  equations  derived  immediately 
implies  that  the  solution  of  the  associated  Riccati  equation  is  strongly  differentiable 
in  time,  and  it  allows  us  to  define  a  "strong”  solution  of  the  Riccati  equation.  Another 
important  consequence  of  this  is  that  the  optimal  control  for  the  linear  quadratic 
regulator  (LQR)  problem  is  continuously  differentiable  if  the  initial  datum  is  sufficiently 
smooth. 

The  Chandrasekhar  equations  for  infinite  dimensional  systems  have  been  discussed 
in  [4]  and  [7]  using  a  Lions-tvpe  framework  [17],  However,  the  equations  derived  in 
[4]  and  [7]  are  satisfied  in  the  distributional  sense.  In  [22],  Sorine  derived  a  set  of 
Chandrasekhar  equations  satisfied  in  a  strong  sense  for  parabolic  systems.  Sorines 
derivation  relied  on  the  analvticitv  of  the  semigroup  and  thus  does  not  apply  to  general 
systems.  Our  approach  differs  from  those  above  in  that  it  uses  an  approximation 
technique.  A  sequence  of  approximating  optimal  control  problems  is  chosen  for  which 
the  Chandrasekhar  equations  may  be  derived  as  in  the  finite  dimensional  case  (see 
[6],  [14],  and  [16]).  Convergence  is  then  established  and  the  appropriate  equations 
are  shown  to  be  satisfied.  In  this  paper,  our  considerations  are  restricted  to  the  LQR 
problem,  but  the  results  are  also  applicable  to  the  Kalman  filtering  problem  [8]. 
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The  spectrum  of  A  is  only  poini  spectrum  and  X  € 
oi>t)  if  and  only  if  del  A(X)  «  0,  where 


A  new  spline-based  scheme  is  developed  for  linear 
retarded  functional  differential  equations  within  the 
framework  of  the  semigroups  on  the  Hilbert  space  R“  a  L3. 
The  approximating  semigroups  preserve  (uniformly)  the 
logarithmic  sectorial  property  (i.e.,  the  differentiability)  of 
the  solution  semigroup.  We  prove  the  convergence  of  the 
scheme  both  in  Wl<*,  R'  *  LJ,  and  using  the  uniform 
differentiability  we  are  able  to  establish  error  estimates 
and  obtain  the  quadratic  convergence  for  a  class  of  initial 
data.  We  also  apply  the  scheme  for  computing  the 
feedback  solutions  to  the  linear  quadratic  regulator 
problem.  The  uniform  differentiability  is  again  the 
essential  basis  for  obtaining  a  complete  theory  for 
convcr&tncc_ot  the  Riccati  solutions. 


0(X)  -  XI  -  L(ck  I).  X  €  C 

We  shall  also  use  the  state-space  X  -  H3(-r,0;P")  for 
system  (1.1).  X  endowed  with  the  inner  product 

<C.tk>x  -  <«0),lK0)>pn  *  <<.tlr>L5 

is  metricaloly  isomorphis  to  dom  A  endowed  with  the 
graph  norm,  which  is  a  Hilberi  space,  because  A  is  closed. 
The  isomorphism  dom  A  -  X  is  given  by 


t(#0).d)  -  «,  (t’(O).e’)  «  dom  A, 
f3<>  -  («0).d).  C  €  X. 


(1.5) 


In  this  paper  we  consider  the  linear  hereditary 
control  system 

x(t)  -  I  A;x(t  ♦  e()  +  f  A(0)x(t  ♦  6)d6  +  Bu(l).  t  >  0, 

i»0  '  -r 

x(0)  -  n,  x(6)  -  d(6)  a.e.  on  (-r.0),  (1.1) 


y(t)  -  Cx(t).  t  »  0, 


where  -r  -  6,  eo  -  0.  x(t)  €  P".  u(t)  €  P“  and  y(t) 

€  P?.  Furthermore,  A(-)  is  an  n»n-matrix  valued  square 
integrable  function  on  J-r.OJ.  It  is  well-known  that,  for 
(h,d)  €  Z  •  P"  *  L,(-r,0;Rn)  and  u  6  L^fO.^R®),  system 

(1.1)  admits  a  unique  solution  x  €  LJ(-i,T;R")  ft  HI(0,T;Rn) 
for  any  T  >  0.  We  define  the  operators  S(t),  t  >  0,  by 

S(t)(rt.d)  -  (x(t),x,),  t  >  0.  (h.d)  €  Z, 

where  x(t)  is  the  solution  of  (1.1)  with  initial  data  (n.d) 
and  u(t)  e  0.  x,  is  the  function  defined  by  x,(8)  • 
x(t  *  6),  -r  <  e  (  0.  S(-)  is  a  Cc-semigroup  on  Z  with 
inifinitesimal  generator  A  given  by  (see  [!],  for  instance) 

dom  A  -  ((h.d)  €  Z|d  €  HV-r.O;^),  n  -  d<0)), 

A(f( 0),d)  -  (L(d).O)  for  («0),d)  €  dom  A. 
where  I  A> 

Ld  -  r  Ajtffj)  *  I  A(6)4<e)de,  t  €  C(-r,0;R"). 

I=C  J  -T 

If  we  define  the  input  operator  B.  FT  -  Z  and  the  output 
operator  C:  Z  -  P®  by 


In  a  recent  paper  [7)  1.  Lasiecka  and  A.  Manitius 
gave  for  the  first  time  optima)  rates  for  the  AV-scheme 
considered  in  (1).  These  estimates  are  essentially  based  on 
uniform  (with  respect  to  the  approximation  parameter) 
differentiability  of  the  approximation  semigroups,  which 
means  that  the  characterization  of  differentiable  semigroups 
in  [E;  Theorem  4.7)  is  uniformly  valid  for  the 
approximating  semigroups. 

In  this  paper  we  devclope  a  scheme  using  first  order 
splines  which  essentially  has  all  the  good  properties  of  the 
AV-scheme.  In  addition  we  are  able  to  establish  error 
estimates  analogous  to  those  in  [7],  As  expected  we  obtain 
quadratic  convergence  for  sufficiently  smooth  data. 
Naturally  uniform  differentiability  is  also  the  essential 
basis  for  our  approach  which  as  far  as  convergence  rates 
are  concerned  is  motivated  by  [7], 


Let  tj*’  -  -k  Jj,  k  -  0, -  •  -  ,N,  and  t*,  »  0.  tji+,  «  -r 

for  N  ■  1,2,- ••  bJ“\  k  -  0,.-  -,N,  denotes  the  usual  first 
order  B-splincs  on  the  interval  [-r,0]  corresponding  to  the 
mesh  t*,  ..,t* 

’**•£.>  f03  Ci  «  e  <  «?• 

B*(e)  -  jKtJ'.j-B)  for  t*  <  e  «  t*,. 

.0  clsewnere. 

Furthermore,  we  put 

Ex  “  x|t*i*  j)-  k  “  1.-  • 


flu  -  (Bu,0),  u  €  P™, 
C(h.d)  -  CD.  (rt,d>)  €  Z. 


(1-2) 


then  (1.1)  is  equivalent  to  the  following  abstract  system  in 
Z: 

z(t)  -  Az(t)  t  £u(t),  t  *  0. 


K0)  -  (n.d),  (1.3) 

y(t)  -  Cz(t),  t  >  0, 

i.e.  a  function  x:  (0.-)  -  P"  is  a  solution  of  (1.1)  if  and 
only  ir  the  function  z(t)  -  (x(t),x,),  t  >  0,  is  a  mild 
solution  of  (1.3)  which  means 

z(t)  «  S(t)( n.d)  ♦  J  S(t-s)£u(s)ds,  t  >  0.  (1.4) 


A  A 

E0n-  (1,0),  E^'  -  (O.E^'l),  k  -  1,  ■  •  -,N, 

where  1  denotes  the  n»n  identity  matrix.  The  following 
spaces  will  be  used  in  the  sequel:  t 

WN  -  span(E*I,  •  •  -.Ejjl)  C  L*(-i.0;R").  • 

A  A 

ZK  -  lev«WK  -  span  (E^,--.,^)  C  Z, 

XK  -  span(Bpl, •  •  -.Ejjl)  C  X, 

ZK  „  flXK  C  dom  A 

It  is  convenient  to  introduce  the  ’basis’  matrices 
E*  -  (E*I.-  •  ,E}jl), 

En  -  (l£,- •  •.!£). 
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ABSTRACT 

It  is  shown  that  the  adjoints  of  a  spline  based  approximation  scheme  for 
delay  equations  do  not  converge  strongly. 
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